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DETECTION  PERFORMANCE  OF  OR-ING  DEVICE  WITH  PRE-  AND 
POST-AVERAGING:  PART  III  -  DETERMINISTIC  SIGNAL 

INTRODUCTION 

BACKGROUND 

The  need  to  process  and  condense  large  amounts  of  data  is 
encountered  frequently  in  modern  Navy  systems  that  employ 
multiple  beams,  frequency  bins,  range  cells,  et  cetera.  One  way 
of  accomplishing  this  goal  is  by  or-ing  a  number  of  inputs  into  a 
single  output  -  that  is,  allowing  only  the  largest  of  a  set  of 
quantities  to  pass  on  for  further  processing  and  completely 
rejecting  the  remainder.  However,  since  this  or-ing  operation  is 
highly  nonlinear,  destroys  information,  and  tends  to  cause  small- 
signal  suppression,  some  pre-averaging  of  the  inputs  to  the 
device  is  often  employed  in  an  effort  to  build  up  the 
signal-to-noise  ratio  (SNR)  prior  to  the  maximum  comparison. 
Furthermore,  at  the  or-ing  output,  some  additional  post-averaging 
used,  this  time  in  an  effort  to  build  up  the  SNR 
before  a  threshold  comparison  is  made  for  the  purpose  of 
declaring  a  signal  present  or  absent  (hypothesis  or  Hq, 

respectively).  The  pertinent  block  diagram  is  displayed  in 
figure  1. 

There  are  N  channels  of  real  input  data  available  for 
processing,  namely,  {x^(t)}  for  n=l:N,  where  time  t  has  been 
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PRE¬ 
INPUTS  AVERAGERS 


POST- 

OR-ING  AVERAGER  OUTPUT 


Figure  1.  Or-ing  with  Pre-  and  Post-Averaging 


normalized  so  that  time-sampling  instant  t  is  an  integer.  Under 
hypothesis  Hq ,  there  is  only  Gaussian  noise  in  all  the  inputs, 
whereas  under  hypothesis  a  signal  is  also  present  in  one 
(unknown)  channel.  The  goal  of  the  processor  in  figure  1  is  to 
determine  signal  presence  with  a  high  detection  probability  P^, 
while  realizing  a  specified  acceptable  low  false  alarm 
probability  P^. 

Each  pre-averager  accumulates  K  statistically  independent 
consecutive  time  samples  of  its  corresponding  input  x^(t), 
yielding  output  >  which  is  then  subjected  to  or-ing  amongst 

N  competitors.  The  or-ing  output  at  time  t  is 

v( t)  =  max{y^( t) , • • • ,y^(t) }  .  (1) 

Finally,  the  post-averager  accumulates  M  samples  of  its  input 
v(t)  and  compares  its  output  w(t)  with  a  fixed  threshold.  It  is 
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presumed  that  the  post-averager  waits  for  all  the  K  input  data 
samples  to  be  summed  and  or-ed  before  an  input  is  received;  then 
another  block  of  K  input  samples  is  summed  and  or-ed,  leading  to 
the  next  input  to  the  post-averager.  That  is,  block  processing 
is  presumed. 

SCOPE 

This  technical  report  constitutes  part  III  of  a  series  of 
three  investigations  of  or-ing  with  pre-  and  post-averaging.  In 
particular,  three  different  input  signal  models  in  additive 
Gaussian  noise  are  of  interest: 

(I)  Random  (Gaussian)  signal, 

(II)  Phase-incoherent  signal,  and 

(III)  Coherent  (deterministic)  signal. 

Results  for  part  I,  the  Gaussian  signal,  have  been  documented  in 
reference  1.  Specifically,  the  details  on  block  and  overlap 
processing,  as  well  as  a  scaling  property  and  the  standard  of 
comparison  for  the  or-ing  system  in  figure  1,  are  presented  here 
along  with  references  to  related  past  results.  Part  II,  the 
phase-incoherent  signal,  was  documented  in  reference  2.  The 
reader  should  be  familiar  with  the  contents  of  both  of  these 
earlier  documents-. 


The  current  report  addresses  part  III,  the  deterministic 
signal  in  noise,  and  presents  numerous  receiver  operating 
characteristics  (ROCs)  that  completely  quaritify  the  performance 
of  the  or-ing  system  in  figure  1.  For  all  three  signal  models 
above,  the  optimum  processors  have  been  determined,  as  shown  in 
appendix  A  of  reference  1.  Thus,  the  exact  losses  that  the 
or-ing  procedure  causes  can  be  accurately  quantified. 

For  the  block  processing  in  figure  1,  expressions  are  derived 
that  allow  for  accurate  evaluation  of  the  false  alarm  probability 
Pj  and  the  detection  probability  for  the  decision  variable 
w(t).  Furthermore,  these  evaluations  are  accomplished  for 
arbitrary  amounts  of  pre-averaging  K,  arbitrary  amounts  of  or-ing 
N,  arbitrary  amounts  of  post-averaging  M,  and  arbitrary  input 
SNRs.  No  approximations  are  involved,  the  analysis  is  not 
limited  to  mean  and  variance  calculations,  and  no  appeal  is  made 
to  the  central  limit  theorem.  Rather,  the  approach  employs  a 
judicious  combination  of  analysis  and  computer-aided  numerical 
calculations.  The  accuracy  of  the  end  result  is  limited  only  by 
the  accuracy  of  the  computer.  An  entire  ROC  can  be  generated  in 
minutes. 
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INPUT  DATA  MODELS 


In  this  section,  gj^  is  a  zero-mean,  unit-variance  Gaussian 
random  variable  (RV). 

DETERMINISTIC  SIGNAL  IN  GAUSSIAN  NOISE 

Here,  in  part  III,  the  signal  channel  input  at  time  k  under 
is  amplitude  variate 

*k  '  9k  ’*k  '  <2> 

where  dj^  is  nonrandom.  The  input  mean  is  3^  =  dj^,  while  the 

input  variance  is  var(Xj^)  =  1.  Thus,  dj^  is  an  input  amplitude 

2  2 

SNR  measure,  while  dj^  (or  dj^/2)  is  an  input  power  SNR  measure. 

For  the  noise  channels,  Xj^  =  gj^,  with  zero  mean  and  unit 
variance.  Thus,  dj^  can  be  interpreted  as  a  per-sample  difference 
in  mean  inputs  that  is  divided  by  the  common  input  standard  devi¬ 
ation.  This  quantity,  dj^,  will  be  called  the  input  deflection. 

When  dj^  =  d  for  all  k  (that  is,  a  common  input  deflection  d 
over  time,  or  a  constant  signal  amplitude  with  time),  the  pre¬ 
averager  output  under  for  the  signal  channel  is 

K  K 

y“n]xk  =  EI!gk  +  Kd  =  y°  +  Kd,  (3) 

k«l  k=l 

where  y°  is  the  corresponding  noise-only  pre— averager  output. 

The  ROCS  can  be  parameterized  by  input  deflection  d  or  input 
power  SNR  measure  d^  (or  d^/2).  More  generally,  the  single 
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parameter  value  =  Ej^  dj^  characterizes  the  performance  of  the 
processor  when  the  individual  input  deflections  {dj^}  are  unequal 
and  arbitrary.  However,  this  observation  is  of  limited  utility 
unless  all  of  the  M  sums  of  one  for  each  post— averaging 

interval,  are  identical. 


OPTIMUM  PROCESSING  OF  AVAILABLE  DATA 


The  optimum  processor  for  signal  model  (2)  has  been  derived 
in  appendix  A  of  reference  1  for  the  case  where  channel  occupancy 
identification  is  required,  as  well  as  for  the  alternative  case 
where  channel  identification  is  irrelevant.  For  all  situations, 
the  optimum  processor  is  essentially  given  by  the  threshold 
comparison 


max 

n=l:N 


Xn(t) 


(4) 


see  equations  (A-13),  (A-19),  and  (A-26)  of  reference  1.  This 
operation  is  equivalent  to  performing  all  pre-averaging  and  no 
post-averaging  in  figure  1,  that  is,  taking  K  =  T,  the  total 
time-bandwidth  product  available,  and  M  =  1.  Doing  so  defers  the 
nonlinearity  (or-ing)  until  after  all  possible  linear  combining 
(pre-averaging)  has  been  accomplished.  This  case  will  be 
thoroughly  evaluated  numerically  and  will  serve  as  the  basis  of 
comparison  for  the  various  combinations  of  pre-averaging  and 
post-averaging,  that  is,  K  <  T  or  M  >  1. 
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PERFORMANCE  ANALYSIS  FOR  BLOCK  PROCESSING 


Due  to  the  pre-averaging,  each  RV  in  figure  1  is  a  sum 

of  K  independent  RVs  with  identical  statistics.  Also,  all  N 
channel  inputs  are  statistically  independent  of  each  other. 

For  or-ing  to  be  present,  N  >  2  is  required  for  the  following 
calculations  and  considerations. 

For  a  noise-only  channel,  let  Pq  be  the  probability  density 
function  (PDF)  of  pre-averager  output  yjj(t),  and  let  Cq  be  the 
corresponding  cumulative  distribution  function  (CDF),  namely, 
Cq(u)  =  Prob{y^(t)  <  u|Hq}.  For  a  signal  present  in  channel  1, 
say,  let  p^  and  c^  be  the  corresponding  PDF  and  CDF  of  y]^(t). 

The  or-ing  output  v(t)  in  figure  1  is 

v(t)  =  max{y^(  t)  ,  .  .  .  ,yj^(  t) }  .  (5) 

Its  CDF  for  signal  present,  that  is,  hypothesis  is 

c^l(u)  =  c^(u)  [Cq(u)]^'^  ,  (6) 

which  leads  to  the  corresponding  PDF  of  v(t)  under  as 
p^]^(u)  =  c^^(u)  =  Pj^(u)  [Cq(u)]^“^  +  Cj^(u)  (N-1)  [Cq(u)]^”^ 

=  [Cg(u)]^"^  [P3^(u)  Cq(u)  +  (N-1)  Pq(u)  Cj^(u)]  .  (7) 

The  signal-absent  PDF  of  v(t)  under  Hq  follows  immediately  as 

Pvo(u)  "  Pq^^^  [Cq(u)]^“^  .  (8) 
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The  corresponding  exceedance  distribution  functions  (EDFS) 
(namely,  e^(u)  =  1  -  c^(u)  for  or-ing  output  v(t))  can  be 
expressed  in  terms  of  the  following  series  expansions,  which  are 
useful  for  very  small  EDF  values: 

N-1 

evi(u)  =  e^(u)  -  c^(u)  H]  [  n  J 

evo'“>  '  -  E  G!)  I-e„(u)l"  .  (10) 

n=l 

The  characteristic  function  (CF)  of  or-ing  output  v(t)  is 
given  by  the  Fourier  transform 

f^(5)  =  exp( i^v( t ) )  =  J  du  exp(i^u)  p^(u)  ,  (11) 

where  either  relevant  form,  p^^^  or  p^Q  above,  is  to  be  used. 
Finally,  for  the  block  processing  in  figure  1,  there  is  no 
overlap  of  the  data  used  in  the  post-averager,  meaning  that  the 
processor  output  w(t)  is  the  sum  of  M  independent  identically 
distributed  RVs.  Therefore,  the  CF  of  the  decision  variable  w(t) 
in  the  figure  is  simply 

fJK)  »  [f^(U]”  .  (12) 

One  final  Fourier  transform  is  required  to  manipulate  CF 
f^(5)  into  the  desired  EDF  e^(u)  of  processor  output  w(t);  see 
references  3  and  4.  This  EDF  will  be  the  false  alarm  probability 
or  the  detection  probability  of  the  or-ing  processor, 
depending  on  whether  PDF  p^Q  or  p^^^  is  used,  respectively. 
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The  fundamental  Fourier  transform  in  equation  (11)  will 
generally  have  to  be  accomplished  numerically  by  means  of  a  fast 
Fourier  transform  (FFT).  If  pre-averager  output  statistics  Pq 
and  Cq  can  be  found  in  closed  form  or  a  readily  computed  form,  a 
cascaded  FFT  approach  will  lead  to  exact  false  alarm  probability 
results  for  the  or-ing  processor  in  figure  1.  If  pre-averager 
output  probability  functions  p^  and  c^  can  also  be  readily 
evaluated,  accurate  detection  probability  results  can  be 
similarly  calculated  as  well.  Some  important  numerical 
considerations  for  this  cascaded  FFT  approach  are  presented  in 
appendix  B  of  reference  1. 

For  the  special  case  of  M  =  1  (no  post-averaging),  processor 
output  RV  w(t)  =  v(t)  in  figure  1,  and  the  corresponding  EDFs 
follow,  for  any  K  and  N,  as 

e^^(u)  =  1  -  c^(u)  tcQ(u)]^"^  ,  (13) 

e^Q(u)  =  1  -  tCQ(u)]”  .  (14) 

Now,  return  to  arbitrary  values  of  K,  N,  and  M,  and 
specialize  the  general  results  above  to  the  deterministic  signal 
model  of  interest  here,  namely,  equations  (2)  and  (3). 
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DETERMINISTIC  SIGNAL 


Define  d  =  /K  d  and  auxiliary  Gaussian  functions 


<|)(x)  =  exp(-x^/2) 


$(x) 


•Ku)  . 


(15) 


Then,  for  constant  signal  amplitudes,  the  requisite  probability 
functions  for  pre-averager  outputs  {yj^(t)}  for  any  K  and  u  are 


Po(u) 

Pl(u) 

Cfllu) 


(2  UK) 


hi?  Jk)  -  tI  ♦(t#)  ' 


( 2  UK ) 


exp 


(- 


(u  -  Kd) 
2K 


Ktk)  '  =  ^(tk  ■  ^] 


(16) 


(17) 


(18) 


Using  equations  (8)  and  (7),  the  PDFs  for  or-ing  output 
v(t)  for  any  K  and  N,  under  Hq  and  are,  respectively, 

PvO^^^  =  A  *(7k)  ^(tk)  '  <19) 

Pvi^^^  =  A  Ktk]  b[A  -  Ktk]  (N-1)  ♦(tI)  *(7^  -  d]]  . 

(20) 

The  corresponding  CFs  are  given  in  integral  form  as 

^vo(^)  =  N  J  dx  exp(i^K^x)  i>(x)  §(x)^~^  ,  (21) 

fvi<^)  =  Jdx  exp(i^K*^x)  §(x)^”^  [<|>(x-d)  $(x)  +  (N-l)  ij)(x)  $(x-d)]. 

(22) 
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Equation  (12)  can  now  be  utilized  to  obtain  the  CFs  of  output 
decision  variable  w(t)  in  figure  1  under  both  hypotheses,  Hg  and 


A  SHORTCUT  IN  TABULATED  RESULTS 

For  a  specified  and  P^,  let  d(K,N,M)  be  the  input 
deflection  required  to  achieve  that  operating  point.  In 
particular,  let 

d^(N,M)  =  d(l,N,M)  (23) 

be  the  input  deflection  required  for  K  =  1.  Also,  observe  from 
the  forms  of  equations  (16)  through  (18)  that  the  effect  of 
adding  K  independent  Gaussian  RVs  is  simply  to  scale  the  input 
deflection  measure  by  /K,  which  means  that  the  required  input 
deflection  satisfies  the  relation 


d(K,N,M) 


d(l,N,M) 


d3^(N,M) 

7k 


(24) 


Therefore,  it  is  only  necessary  to  evaluate  ROCs  and  tabulate 
numerical  results  for  the  case  of  K  =  1;  N  and  M  can  be  general 
ihtegers.  This  observation  significantly  reduces  the  amount  of 
computations  and  graphs  that  need  to  be  investigated  and 
catalogued.  In  the  following,  the  ROCs  will  be  computed  solely 
for  K  =  1  and  labeled  with  the  corresponding  values  of  input 
deflection  dj^(N,M).  However,  these  labels  could  be  replaced  by 
/K  d(K,N,M)  for  any  K,  and  the  same  ROCS  would  still  apply  for 
any  selected  value  of  integer  K. 
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SPECIAL  CASES 


When  M  =  1  (no  post-averaging),  the  EDFs  in  equations  (13) 
and  (14)  can  be  combined  with  the  results  in  equation  (18)  to 
immediately  yield  closed-form  expressions  for  and  for  any 
K,  N,  and  input  deflection  d,  namely, 

Pj  =  1  -  [f(u//K)]^  ,  P^  =  1  -  f(u//K-  d)  [§(u//K)]^"^  .  (25) 


Here,  d  =  d  /K.  Furthermore,  these  equations  can  be  solved 
explicitly  for  the  required  threshold  u  and  the  required  input 
deflection  d  to  realize  a  specified  P^  and  P^: 


u  =  /K  1(a) 


d(K,N,l) 


5(a)  -  5(ae)  di(N,l) 
7k - —7k - 


(26) 


where  $  is  the  inverse  i  function,  and  the  auxiliary  parameters 
a  and  p  are  defined  as 


a  =  a(N)  s  (1  -  P^)^^^  , 


(27) 


For  N  =  1,  no  or-ing, 
results  above  simplify  to 


and  arbitrary  K  and  M,  under  the 

system  output  PDF 


Vl 


(u)  = 


(2iiT)’ 


exp 


(u  -  Td) 
2T 


■)  - 


T  =  K  M 


(28) 


and  output  EDF 


Pa(v) 


.  I  du  p„^(u)  .  ♦(qd  -  5)  -  g 


=  =  (KM)*^ 


(29) 
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Obviously,  for  N  =  1,  only  the  product  KM  matters,  at  least  for 
unity  weights  in  both  the  pre-  and  post-averagers. 

For  the  special  case  of  N  =  2,  the  integrals  in  equations 
(21)  and  (22)  can  be  evaluated  in  terms  of  the  error  function  of 
complex  argument  w(  );  see  reference  5,  chapter  7.  The  CFs 
follow  as 

f^Q(5)  =  exp(-Y^)  w(y)  (30) 

and 

1  2  2 

fvi(^)  =  2  +i2Yv)  [w(Y+iv)  +  w(Y-iv)]  ,  (31) 

k  k 

where  y  =  K^^/2  and  v  =  K^d/2.  The  amount  of  pre-averaging  K  is 
arbitrary.  The  CFs  of  decision  variable  w(t)  are  given  by  the 
M-th  powers  of  the  above  two  results  under  Hq  and  H^, 
respectively,  when  N  =  2. 


NUMERICAL  APPROACH 

When  the  CF  f^(^)  of  or-ing  output  v(t)  can  be  derived  in 
closed  form,  such  as  in  the  special  cases  (30)  and  (31),  f^(^) 
can  be  used  directly  in  equation  (12)  to  find  the  closed-form  CF 
of  decision  variable  w(t).  This  route  is  preferred  because  CF 
f^(E)  can  decay  rather  slowly  with  5  for  some  values  of  K  and  N, 
making  the  use  of  equation  (11)  difficult. 
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Therefore,  closed-form  results  for  CF  f^(^),  when  available, 
can  be  very  useful  in  terms  of  avoiding  some  numerical  problems 
in  the  discrete  Fourier  transform  of  p^(u)  according  to  equation 
(11).  When  this  is  not  possible,  the  closed— form  results  for  the 
or-ing  output  PDFs  given  in  equations  (19)  and  (20)  are  used, 
instead,  in  the  cascaded  FFT  procedure  outlined  in  equations  (11) 
and  (12). 


14 


QUANTITATIVE  PERFORMANCE  RESULTS 


REQUIRED  INPUT  DEFLECTIONS 

This  section  gives  specific  quantitative  detectability 
results  with  various  amounts  of  or-ing  N  and  post-averaging  M, 
while  the  amount  of  pre-averaging  is  kept  at  K  =  1.  In 
particular,  54  ROCs  are  presented  in  appendix  A  for  the  cases  of 

N  =  1,  2,  4,  8,  16,  32, 

M  =  1,  2,  4,  8,  16,  32,  64,  128,  256. 

A  standard  operating  point  (SOP)  is  defined  here  as 
Pj  =  lE-3  and  P^  =  0.5,  while  a  high-quality  operating  point 

(HOP)  is  defined  as  P^  =  lE-6  and  P^  =  0.9.  The  required  input 

deflection  values  dj^(N,M)  for  the  SOP  and  HOP  have  been  extracted 
from  the  ROCs  in  appendix  A  and  are  presented  in  tables  1  and  2, 
respectively.  These  results  for  K  =  1  can  be  used  to  find  the 
required  input  deflection  for  general  integer  K  according  to 

d^ (N,M) 

d(K,N,M)  =  -  .  (32) 


In  conjunction  with  equation  (32),  the  input  deflections  in 
tables  1  and  2  form  the  basis  for  the  plots  in  figures  2  through 
9  for  both  the  SOP  and  HOP  cases  of  KM  =  4,  16,  64,  and  256. 
Figures  10  and  11  are  identical  to  figures  8  and  9,  except  that 
the  ordinates  are  now  in  decibels  according  to 
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Table  1. 

Required 

Input 

Deflection 

d,(N,M)  for 

=  lE-3 

,  -  0. 

5,  K  = 

1,  Deterministic  Signal 

N 

2 

4 

8 

16 

32 

M 

1 

3.09 

3.29 

3.48 

3.66 

3.83 

4.00 

2 

2.19 

2.44 

2.67 

2.90 

3.10 

3.30 

4 

1.55 

1.82 

2.08 

2.33 

2.56 

2.77 

8 

1.09 

1.37 

1.63 

1.89 

2.13 

2.35 

16 

0.773 

1.02 

1.28 

1.53 

1.78 

2.01 

32 

0.546 

0.756 

0.994 

1.24 

1.49 

1.72 

64 

0.386 

0.557 

0.766 

0.998 

1.24 

1.47 

128 

0.273 

0.407 

0.584 

0.794 

1.02 

1.25 

256 

0.193 

0.296 

0.440 

0.623 

0.833 

1.05 

Table  2.  Required  Input  Deflection  d^(N,M)  for 
Pf  =  lE-6,  =  0.9,  K  =  1,  Deterministic  Signal 


M 


N 


16 


32 


1 

2 

4 

8 

16 

32 

64 

128 

256 


6.03 

6.17 

6.31 

6.44 

6.57 

6.69 

4.27 

4.46 

4.64 

4.82 

4.98 

5.15 

3.02 

3.27 

3.50 

3.71 

3.91 

4.11 

2.13 

2.42 

2.68 

2.93 

3.15 

3.37 

1.51 

1.80 

2.08 

2.34 

2.58 

2.81 

1.07 

1.35 

1.62 

1.89 

2.14 

2.37 

0.754 

1.00 

1.26 

1.53 

1.78 

2.02 

0.533 

0.742 

0.980 

1.23 

1.48 

1.72 

0.377 

0.546 

0.753 

0.986 

1.23 

1.46 

dB  =  20  log(d),  since  d  is  an  amplitude  SNR  measure.  Additional 
figures  could  also  be  plotted  for  the  intermediate  cases  of 
KM  =  2,  8,  32,  and  128,  if  desired,  by  use  of  tables  1  and  2. 
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Figure  2.  Required  Iijput  Deflecticn  d(K,N,M) 

Pj  =  lE-3,  =  0.5,  KM  =  4,  Deterministic  Signal 


Figure  3.  Required  Input  Deflection  d(K,N,H) 

Pj  =  lE-6,  Pj  =  0.9,  KM  =  4,  Deterministic  Signal 


Figure  4.  Required  Input  Deflection  d(K,N,M) 

Pj  =  lE-3,  =  0.5,  KM  =  16,  Deterministic  signal 


Figure  6.  Required  Input  Deflection  d(K,N,H) 

Pf  =  lE-3,  =  0.5,  KM  =  64,  Deterministic  Signal 


1.1 


0.9 


•D 

C 

,o 


"S 

Q 


0.7 


0.5 


0.3 


0.1 


M 


Figure  8.  Required  Input  Deflection  d(K,N,M) 
for  =  lE-3,  =  0.5,  KM  =  256,  Deterministic  Signal 


Figure  9.  Required  Input  Deflection  d(K,N,M) 
for  Pj  =  lE-6,  P^  =  0.9,  KM  s  256,  Deterministic  Signal 


20 


1  2  4  8  16  32  64  128  256 

M 

Figure  10.  Required  Ii^t  Deflection  d(K,N,M)  in  Decibels 
for  Pj  =  lE-3,  =  0.5,  KM  =  256,  Deterministic  Signal 


1  2  4  8  16  32  64  128  256 


M 


Figure  11.  Required  Input  Deflection  d(K,N,M)  in  Decibels 
for  Pj  =  lE-6,  Pj  =  0.9,  KM  =  256,  Deterministic  Signal 


The  accurate  deflection  values  in  tables  1  and  2  were 
interpolated  from  the  ROCs  in  appendix  A  while  the  false  alarm 
and  detection  probability  numbers,  and  P^,  were  still  resident 
in  the  computer;  eyeball  interpolation  from  the  plotted  ROCs 
cannot  be  accomplished  this  accurately. 

It  can  be  seen  from  figures  2  through  9,  as  expected,  that 
the  required  input  deflection  d  increases  monotonically  with  the 
number  (N)  of  channels  or-ed  if  K  and  M  are  held  fixed.  Also, 
the  required  input  deflection  increases  monotonically  with  the 
amount  (M)  of  post-averaging  employed  if  N  and  product  KM  are 
held  fixed;  alternatively,  the  required  input  deflection 
decreases  with  the  amount  (K)  of  pre-averaging  employed  if  N  and 
product  KM  are  held  fixed.  The  exact  rates  can  be  determined 
from  figures  2  through  9  and  the  ROCs  in  appendix  A. 


TRADEOFF  BETWEEN  PRE-  AND  POST-AVERAGING 

For  a  fixed  amount  of  total  time-bandwidth  product,  that  is, 

T  =  KM  =  constant,  it  is  possible  to  shift  some  of  the  pre¬ 
averaging  to  post-averaging,  or  vice  versa.  For  example,  if 
KM  =  16,  this  shift  could  be  accomplished  by  taking  K  =  16  and 
M  =  1,  which  corresponds  to  no  post-averaging,  and  then  switching 
to  K  =  1  and  M  =  16,  which  corresponds  to  no  pre-averaging. 
Alternatively,  intermediate  values  such  as  K  =  4  and  M  =  4  could 
be  used,  still  keeping  KM  =  16. 
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In  all  cases,  increasing  M  at  the  sake  of  K  requires  larger 
input  deflections  to  maintain  the  same  operating  point;  see 
figures  2  through  9.  This  is  because  the  optimum  processor 
conducts  all  pre-averaqing  before  any  nonlinear  operations;  see 
appendix  A  of  reference  1.  For  example,  to  maintain  the  SOP  when 
KM  =  16,  table  1  and  equation  (32)  indicate  that  the  input 
deflection  must  be  scaled  by  factors  of  1.24,  1.47,  1.67,  1.86, 
and  2.01  for  N  =  2,  4,  8,  16,  and  32,  respectively,  if  the  switch 
is  made  from  all  pre-averaging  to  all  post-averaging. 
Alternatively,  to  maintain  the  HOP,  the  corresponding  scaling 
ratios  from  table  2  and  equation  (32)  are  smaller  (1.17,  1.32, 
1.45,  1.57,  and  1.68  for  N  =  2,  4,  8,  16,  and  32,  respectively). 

For  a  HOP  such  as  P^  =  lE-6  and  P^  =  0.9,  larger  input 
deflections  are  naturally  required  to  achieve  the  higher  level  of 
performance.  At  these  larger  input  deflections,  the  or-ing 
process  is  more  frequently  dominated  by  the  signal-bearing 
channel.  Therefore,  as  the  number  (N)  of  channels  or-ed 
increases,  the  required  increases  in  input  deflection  are  less 
severe  than  for  a  lower  quality  operating  point,  such  as  SOP 
Pj  =  lE-3  and  P^  =  0.5.  For  example,  for  K  =  16  and  M  =  1, 
increasing  N  from  2  to  32  requires  a  gain  factor  of  1.22  in  the 
deflection  to  maintain  the  SOP,  whereas  only  a  gain  factor  of 
1.08  is  needed  at  the  HOP. 
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A  related  observation  is  that,  for  a  given  configuration 
(fixed  K  and  M) ,  the  relative  increase  in  input  deflection 
required  to  maintain  the  SOP  is  more  severe  for  the  larger  M 
values  as  the  amount  of  or-ing  increases.  Thus,  for  KM  =  16,  as 
N  increases  from  2  to  32,  a  deflection  factor  increase  of  1.22 
suffices  to  maintain  the  SOP  for  M  =  1,  whereas  a  factor  of  1.97 
is  required  for  M  =  16.  The  corresponding  relative  increases  at 
the  HOP  are  1.08  for  M  =  1  versus  1.56  for  M  =  16,  as  N  increases 
from  2  to  32.  The  behaviors  are  similar  for  the  other  values  of 
KM,  although  the  required  input  deflection  values  are  smaller  as 
KM  increases  because  of  the  additional  observation  times 
allowed . 

Inspection  of  figures  10  and  11  reveals  that  the  decibel 
or-ing  losses  in  the  input  deflection  (20  log(d))  are  much  larger 
than  for  the  random  signal  treated  in  part  I  (reference  1)  and 
the  phase-incoherent  signal  in  part  II  (reference  2).  The  reason 
is  that  if  KM  is  held  fixed  and  M  increases,  then  K  is 
decreasing.  But,  since  the  pre-averager  is  now  (part  III) 
performing  coherent  (amplitude)  addition  rather  than  incoherent 
addition  (that  is,  sums  of  envelopes  squared  in  parts  I  and  II) , 
any  decrease  in  K  is  felt  more  strongly.  So,  the  tradeoff  from 
pre-averaging  (larger  K)  to  post-averaging  (larger  M,  smaller  K) 
is  more  detrimental  to  performance  for  a  deterministic  signal 
than  for  a  random  signal  or  a  phase-incoherent  signal. 
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SUMMARY 


For  a  deterministic  signal  in  Gaussian  noise,  closed-form 
expressions  have  been  derived  for  the  PDFs  at  the  output  of  a 
pre-averager  followed  by  an  or-ing  operation.  These  forms  have 
been  numerically  subjected  to  a  cascaded  pair  of  Fourier 
transforms  to  yield  accurate  results  for  the  EDFs  at  the  output 
of  a  post-averager,  under  both  hypotheses  Hq  and  for  numerous 
values  of  the  parameters  of  the  complete  processor  (figure  1). 
Such  results  enable  investigation  of  the  or-ing  processor  for 
false  alarm  probabilities  P^  in  the  range  of  lE-6  and  smaller; 
there  is  no  need  to  resort  to  lengthy  simulations. 

Numerous  ROCs  have  been  generated  that  enable  users  or  system 
designers  to  quickly  assess  the  losses  to  be  expected  from 
employing  or-ing  in  their  processors.  Also,  quantitative 
evaluation  of  the  tradeoffs  between  pre-averaging  and  post¬ 
averaging  has  been  conducted.  Finally,  because  the  enclosed 
tabulations  will  undoubtedly  not  cover  all  cases  of  practical 
interest,  a  MATLAB  program  for  the  evaluation  of  additional  ROCs 
is  presented  in  appendix  B. 

The  analytical  approach  utilized  here  can  be  extended  to 
include  quantizers  on  the  inputs  in  figure  1,  and  accurate  ROCs 
can  still  be  obtained  for  decision  variable  w{t).  In  fact,  there 
can  be  arbitrary  amounts  of  quantization  (L  levels),  and  the 
input  PDFs  can  be  completely  arbitrary.  This  procedure  will  be 
documented  in  a  future  report. 
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APPENDIX  A  -  ROCS  FOR  DETERMINISTIC  SIGNAL 


This  appendix  contains  the  ROCs  for  or-ing  with  pre-  and 
post-averaging  and  a  deterministic  input  signal.  In  keeping  with 
the  discussion  in  equations  (23)  and  (24),  the  size  of  the 
pre-averager  is  set  at  K  =  1,  without  loss  of  generality.  The 
parameters  N  and  M  are  varied  over  the  ranges 

N  =  1,  2,  4,  8,  16,  32, 

M  =  1,  2,  4,  8,  16,  32,  64,  128,  256, 

giving  a  total  of  54  ROCs.  The  results  in  tables  1  and  2  for  the 
SOP  and  the  HOP  were  extracted  from  these  ROCs.  Results  for 
other  operating  points  can  also  be  obtained  from  these  ROCs,  if 
desired. 

In  the  following  figures,  the  ROCs  have  been  computed  solely 
for  K  =  1  and  labeled  with  the  corresponding  values  of  input 
deflection  dj^(N,M).  However,  these  labels  could  be  replaced  by 
/K  d(K,N,M)  for  any  K,  and  the  same  ROCs  would  still  apply  for 
any  selected  value  of  integer  K.  See  the  discussion  on  page  11. 

To  maximize  the  clear  region  of  the  aliased  approximations  to 
the  CFs  and  EDFs  obtained  by  FFTs  (according  to  references  3  and 
4),  a  procedure  for  using  shifted  RVs  is  presented  in  appendix  C. 
This  alternative  is  very  useful  when  a  wide  range  of  different 
input  SNR  levels  is  of  interest. 
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Figure  A-1.  RDCs  for  K=l,  N=l,  M  =  1  Figure  A-2.  POCs  for  K 


Figure  A-3.  ROCS  for  K=l,  N=l,  M  =  4  Figure  A-4.  ROCs  for  K 


Figure  A-5.  FOCs  for  K  =  1,  N  =  1,  M  =  16  Figure  A-6.  POCs  for  K 
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Figure  A-7.  KOCs  for  K  =  1,  N  =  1,  M  =  64  Figure  A-8.  ROCs  for  K  =  1,  N  =  1,  M  =  128 
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n  =  2048 


660 


A- 8 


Figure  A-13.  ROCs  for  K=l,  Ns2,  M=8  Figure  A-14.  POCs  for  K 


Figure  A-15.  POCs  for  K  =  1,  N  =  2,  M  =  32  Figure  A-16.  ROCs  for  K 
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for  K  =  1,  N  =  2,  M  =  128  Figure  A-18.  ROCs  for  K  =  1,  N  =  2,  M  =  256 


Figure  A-19.  POCs  for  Ksl,  N  =  4,  Figuro  A-20.  POCs  for 
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Figure  A-21.  ROCs  for  K  =  l,  N  =  4,  M  =  4  Figure  A-22.  POCs  for  K 
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Figure  A-25.  POCs  for  K  =  1,  N  =  4,  M  =  64  Figure  A-26.  ROCS  for  K  =  1,  N  =  4,  M  =  128 
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Figure  A-29.  POCS  for  K=l,  N=8,  M  =  2  Figure  A-30.  ROCs  for 


Figure  A-31.  FOCS  for  K=l,  N=8,  M  =  8  Figure  A-32.  ROCs  for  K 
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Figure  A-33.  POCs  for  K  b  i,  n  »  8,  H  s  32  Figure  A-34.  ROCs  for  K 
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Figure  A-37.  POCs  for  K  =  1,  N  =  16,  M  =  1  Figure  A-38.  ROCs  for  K 
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Figure  A-41,  ROCS  for  K  =  1,  N  =  16,  H  =  16  Figure  A-42.  ROCs  for  K 
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Figure  A-45.  POCs  for  K  «  1,  N  =  16,  M  =  256  Figure  A-46 
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Figure  A-47.  ROCs  for  K  =  1,  N  =  32,  M  =  2  Figure  A-48.  ROCs  for  K 
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Figure  A-53.  ROCs  for  K  =  1,  N  =  32,  M  =  128  Figure  A-54.  EOCs  for  K  =  1,  N  *  32,  M  =  256 


APPENDIX  B  -  MATLAB  PROGRAM  FOR  EVALUATION  OF 
RECEIVER  OPERATING  CHARACTERISTICS 


Every  figure  in  each  set  of  ROCs  in  appendix  A  has  the  three 
parameter  values  listed  on  it  that  were  used  to  generate  those 
specific  curves.  They  are  the  bias  b,  sampling  increment  du,  and 
FFT  size  n.  Because  an  additive  constant  to  the  decision 
variable  in  figure  1  does  not  change  the  ROC  (reference  1),  bias 
b  was  added  to  RV  w(t)  and  chosen  so  as  to  make  the  sum  variable 
as  small  as  possible,  but  never  negative  (within  roundoff  error). 
This  procedure  maximizes  the  unaliased  region  when  the  final  FFT 
is  conducted  from  the  CF  domain,  equation  (12),  to  the  EDF  domain 
for  system  output  w(t),  as  well  as  eases  the  requirements  on 
sampling  increment  du  and  FFT  size  n. 

Sampling  increment  du  is  applied  in  the  u  domain  to  equations 
(19)  and  (20);  increment  du  must  be  taken  small  enough  so  that 
the  approximation  to  CF  f^(^)  in  equation  (11)  is  not  aliased 
significantly  in  the  ^  domain.  Finally,  FFT  size  n  must  be  taken 
large  enough  so  that  the  u-domain  span,  n  du,  in  the  final  FFT  of 
equation  (12),  from  the  CF  to  the  EDF,  guarantees  insignificant 
aliasing.  Complete  details  on  this  cascade  procedure  are 
presented  in  appendix  B  of  reference  1. 
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All  three  parameter  values  (b,  du,  n)  have  been  chosen  by  a 

trial  and  error  procedure,  in  which  intermediate  CF  is 

observed  in  the  ^  domain  and  output  EOF  e  (u)  is  observed  in  the 

w 

u  domain.  Repeated  trials  led  to  the  tight  parameter  values 
listed  on  each  set  of  ROCs.  These  listed  parameter  values  can 
serve  as  starting  or  reference  points  for  the  evaluation  of 
additional  ROCs  with  different  values  of  K,  N,  M,  and  d  that  are 
of  interest  to  a  user.  The  following  listing  of  the  MATLAB 
program  corresponds  to  figure  A-45. 

clear,  clf  %  Peort  III,  mm  TR  11,248 

b=-380;  %  Additive  bias  to  w 

du=.12;  %  Sarrpling  iacranent  ia  u 

n=2^11;  %  FFT  size 

kc=l;%FIXED  %  K,  amount  of  pre-averaging 
nc=16;  %  N,  amount  of  or-ing 

mc=256;  %  M,  amount  of  post-averaging 

dnin=.5;  %  Miniirtm  deflection  d 

dinc=.05;  %  Incrainent  in  deflection  d 

num=17;  %  Number  of  ROCs 

disp(  'b  du  n;  kc  nc  me:  ' ) 
di£p(  [b  du  n;  kc  nc  me] ) 

xg=[le-6  le-5  le-4  .001  .002  .005  .01  .02... 

.05  .1  .2  .3  .4  .5]; 

Yg=[le-6  le-5  le-4  .001  .002  .005  .01  .02... 

.05  .1  .2  .3  .4  .5  .6  .7  .8  .9  .95  .98  .99]; 
3^=phiinv(xg)  ; 
ygi^EMinv(yg)  ; 

[Xg,Yg]=^teshgrid(xg,Yg)  ; 

]yfegcf=zeros  (n,  1)  ; 

EDF=zeros  (n,  1)  ; 

Pd=zeros(n,nurrH-l)  ; 

El=n-1;  n2=n/2;  n3=n2+l;  ncl=nc-l; 

pn=pi/n;  d2=du/2;  dxi=2*pi/(n*du); 

arg=pin.*  (I:n2)  '  ; 

rsinc=arg.  /sin(arg)  ; 

ex=e>5)  (b*dxi*  (0  :n2 )  '  *i)  ; 

dispC  isnr  edf(O)-l') 
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for  israo=0:ri]Liri 

d=3±nin+diQc*  (isnr-l)  ; 
if(isiiD=0),  d=0;  ead 
X=zeros  (n,  1)  +zeros  (n,  1)  *i; 
meanv^O;  cdfk=0;  area=0; 
u0=phiiiiv'(10^  (-16/nc) ) ;  k;=floor  (uO/du)  ; 
vMle(cdfk  <  .5  I  area  >  le-20) 
k=k+l;  \jte:du*k;  u2=uk+d2; 
cO=iM(ij2)  ; 
if(d=0),  cl=cO; 
else  cl=phi  (u2-d)  ; 
end 

cdfo=cdfk;  cdfk=cl*c0^cl; 
area=cdfk-cdf o ; 
j=mod(k,n)  ; 

X  ( j+1)  =X  ( j+1) +area; 
rneanv^:Ttiear]v-+area*ijk;  %  lik,  not  u2 

end 

X=fftgreen(X)  ; 

Magcf=loglO  (nagsq(X)  +le-50)  *.5; 

plot(]yiagcf , 'r' ) 

axis([l  n+1  -16  0]);  grid  on 

pai:ise(l) 

X  (2 : n3 )  =conj  (X  (2 : n3 ) ) .  *rsinc  ( 1  :n2 )  ; 
X(n2+2:n)=0; 

X=X.'irc; 

X(l:n3)=X(l:n3).*ex; 

itieanw=qiieanv*inc+b; 

X(1)=0; 

X(2:n3)=X(2:n3)  ./(l:n2) 

X=fftgreen(X)  ; 
a=.5-nnearM/  (n*dLi)  ; 
edf  0=a+iinag  (X  ( 1 ) )  /pi  ; 
disp(  [isnr  edf  0-1] ) 
k=(l:n) 

edf=a-(k-l)  ./n+ionag(X(k) )  ./pi; 

X=edf; 

EDF=loglO  (abs  (edf)  +le-30) ; 
plot(BDF,  'r' ) 

axis([l  n+1  -18  0]);  grid  on 
paused) 


edf=real(X(k) ) ; 
edf  ( edf ,  l-le-12 )  ; 

edf=inax(edf ,  le-12)  ; 

Pd( : ,  isnr+1)  =phiin[v(edf )  ; 
and 
be^ 
pause 
clf 

hold  cn 

set(gcf, 'E^perPosition' , [.25  .25  8  10.5] ) 
plot(xg,Yg,  'k') 
plot(Xg,yg,  'k') 

plot(3^,xg,  'k')  %  zero  SNR  ROC 
plot{Pd( :  ,1)  ,Pd( : ,  [2:nurrH-l] ) ,  'k' ) 
axis([xg(l)  xg(14)  ygd)  yg(21)]) 
axis  off 
v\iiLle  1 

T=irpj.t(  'pf  pd  isnr:  ' )  ; 
if(T(l)=0),  break,  end 
pft=phiinv(T(l) )  ; 
pdt=phiinv(T(2) )  ; 
isnr=T(3)  ; 
for  j=0:riL 

if  (Pd(j+l,l)<pft) ,  break,  end 
ori 

xl=Pd(j+l,l) ; 
x2=Pd(j,l); 
as=Pd(j+l, isnr+1) ; 
bs=Pd(j+l,isnr+2) ; 
cs=Ed( j , isnr+1) ; 
ss= (cs-as) / (x2-xl) ; 
al=pdt-ss*pf  t  ; 
fs= (al+ss*xl-as) / (bs-as) ; 
dc=dmin+dinc*  (isnr-l+fs)  ; 
disp ( [T  f s  dc] ) 
eri 

%  function  Y^hiinv(x)  %  0  <  x  <  1 
%  y=1.414213562373095*erfiw(2*x-l)  ; 

%  function  y=phi(x) 

%  . 5*erfcore (-x* . 7071067811865476 , 1) ; 


APPENDIX  C  -  SHIFTED  DECISION  VARIABLES 


Suppose  that  the  ROCs  for  decision  variables  Wq  and  under 
hypotheses  Hq  and  respectively,  are  of  interest.  That  is,  a 
plot  of  exceedance  probability  =  Prob(Wj^  >  kt)  versus 

exceedance  probability  P£(k)  *  Prob(wQ  >  kA)  is  desired.  Let 
shifted  decision  variables 

Zq  =  Wq  +  IqA  ,  Zj^  =  +  I^A  ,  (C-1) 

where  Iq  and  are  integers.  Then, 

Pj(k)  =  Prob(zQ  >  (k  +  Iq)A)  s  QQ(k  +  Iq)  , 

P^(k)  =  Prob(z^  >  (k  +  Ij^)A)  2  ( k  +  I^)  ,  ( C-2 ) 

where  exceedance  probabilities  Qq  and  are  computed  and  stored 
for  decision  variables  Zq  and  under  hypotheses  Hq  and  H^^, 
respectively.  That  is,  probabilities  QQ(n)  =  Prob(zQ  >  nA)  and 
Q^(n)  =  Prob(Zj^  >  nA)  are  available. 

Then,  to  obtain  the  ROC  for  variables  Wq  and  ,  plot 
Q^(k  +  Ij^)  versus  QQ(k  +  Iq),  or,  equivalently,  Qj^(k  +  1^^  -  Iq) 
versus  QQ(k).  The  integers  Iq  and  1^^  should  be  taken  with  as 
large  negative  values  as  possible,  but  still  keeping  the  shifted 
RVs  Zq  and  z^  positive.  This  procedure  enables  maximum  use  of 
the  fundamental  "cleared  region"  in  the  EDF  plot  obtained  by  an 
FFT  procedure;  see  references  3  and  4.  In  fact,  for  best 
results,  can  be  a  function  of  (SNR  measure)  deflection  d. 
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